Uniform Modulus of Continuity of Random Fields 



Yimin Xiao * 
Michigan State University 

June 2, 2009 



Abstract 

A sufficient condition for the uniform modulus of continuity of a random field X = 
{X(t),t € 1*} is provided. The result is applicable to random fields with heavy-tailed 
distribution such as stable random fields. 
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1 Introduction 

Sample path continuity and Holder regularity of stochastic processes have been studied by many 
authors. The celebrated theorems of Kolmorogov [cf. Khoshnevisan (2002)] and Garsia (1972) 
provide uniform modulus of continuity for rather general stochastic processes and random 
fields. For Gaussian processes, a powerful chaining argument leads to much deeper results; see 
the recent books of Talagrand (2006), Marcus and Rosen (2006) and Adler and Taylor (2007). 
Some of the results on uniform modulus of continuity have been extended to non-Gaussian 
processes provided the tail probabilities of the increments of the process decay fast enough [see 
Csaki and Csorgo (1992), Kwapieh and Rosihski (2004)]. 

In recent years, there has been increasing interest in sample path continuity of stochastic 
processes with heavy-tailed distributions such as certain infinitely divisible processes including 
particularly stable processes. See Samorodnitsky and Taqqu (1994), Marcus and Rosihski 
(2005) and the references therein for more information. Since many of such processes do not 
have finite second moment, the aforementioned results on the uniform modulus of continuity 
do not apply. There have been only a few results on uniform modulus of continuity for stable 
processes; see Bernard (1970), Takashima (1989), Kono and Maejima (1991a, 1991b), Bierme 
and Lacaux (2009), Ayache, Roueff and Xiao (2007, 2009). 
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The objective of this paper is to study uniform modulus of continuity of a real-valued 
random field X = {X(t),t G W N } with heavy-tailed distributions. In Section 2, we modify 
the chaining argument to prove a general result on uniform modulus of continuity for X. 
Compared with Kolmogorov's continuity theorem, Dudley's entropy theorem or the theorems 
in Csaki and Csorgo (1992), our result does not assume exponential- type tail probabilities nor 
higher moments, hence it can be applied to wider classes of random fields. In Section 3, we 
introduce the notion of maximal moment index for a sequence of random variables and study 
its basic properties. We show that the maximal moment index is closely related to the uniform 
Holder exponent of the random field X. In Section 4, we apply the results in Sections 2 and 3 
to stable random fields. 

Throughout this paper, we will use K to denote an unspecified positive constant which 
may differ in each occurrence. Some specific constants will be denoted by Ki,K2, — 
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2 A general result on uniform modulus of continuity 

Let X = {X(t),t G T} be a real- valued random field on a compact metric space (T, d). In this 
section we establish a general result on uniform modulus of continuity of X. Similar to the 
proofs of Kolmogorov's continuity theorem and Dudley's entropy theorem [see Khoshnevisan 
(2002), Talagrand (2006), Marcus and Rosen (2006) or Adler and Taylor (2007)], our method 
is also based on the powerful chaining argument. However, in order to be able to apply the 
method to random fields with heavy-tailed distributions such as stable random fields, we will 
make two modifications. One is that we will rely more on the geometry of the parameter set 
T and choose approximating chains more carefully so that, at the n-th step of the chain, there 
are Ko N n pairs of random variables, where N n is the 2 _n -entropy number of T and kq is a 
constant, instead of N n N n _\ pairs as in Talagrand (2006) or Adler and Taylor (2007). This is 
important for deriving correct rate function for the modulus of continuity for processes with 
heavy-tailed distributions. The second modification is, instead of trying to control the tail 
probabilities of the increment of X(t) over every step of the chain as in the aforementioned 
references, we control the 7-th moments of the maximum increments over the steps. Here 
7 G (0, 1] is a constant. This latter argument extends an idea of Kono and Maejima (1991b) 
who studied the modulus of continuity of an //-self-similar a-stable process with stationary 
increments and satisfying a G (1,2) and ^ < H < 1. 

We assume that there is a sequence {D n ,n > 1} of finite subsets of T satisfying the following 
conditions: 

(i). There exists a positive integer kq depending only on (T,d) such that for every r n+ i G 
D n+ i, the number of points r n G D n satisfying d(r n+ i,r n ) < 2~™ is at most kq. Such 
a point r n is called a D„-neighbor of T n+ \. We denote by O n (r n+ i) the set of all D n - 
neighbors of T n+ \. 
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(ii). [Chaining property] For every s,t G T with d(s,t) < 2~ n , there exist two sequences 
{ t p( s ),P — n } an d {r p (t),p > n} such that r n (s) = r n (t) (i.e., the two chains have the 
same starting point) and, for every p > n, both r p (s) and r p (t) are in D p , 

d(r p (t), t) < 2-p, d(r p (s), s) < 2-p 

and 

r p (t) g O p (T p+1 (t)), t p (s) g O p (t p+1 (s)). 

If, in addition, s G D := UfeLi ^fc ( or t £ D), then there exists an integer q > n such 
that Tp(s) = s (or r p (t) = t) for all p > q. 

Condition (ii) implies that, for every n, T is covered by d-balls with radius 2~ n and centers 
in D n ; and the collection of all the points in {D n ,n > 1} is dense in T, i.e., Un=i D n = T. 

When T is a closed interval in M. N and d is a metric which is equivalent to the i°° metric 
\s — t\oo = maxi<j<7v \ sj —tj\, it is convenient to choose a sequence {D n ,n > 1} satisfying the 
above conditions. For example, if T = [0, 1]^ and d is the £°° metric, then, for any integer 
n > 1, one can take D n to be the collection of all vertices of dyadic cubes of order n which are 
contained in T. Note that in this case, {D n ,n > 1} is increasing (i.e., D n C D n+ \ for every 
n > 1) and for every r n G D n , the number of -D n _i-neighbors of r n is at most kq := 3 N — 1. 

The following is the main result of this section. 

Theorem 2.1 Let X = {X(t),t G T} be a real-valued random field on a compact metric space 
(T,d) and let {D n ,n > 1} be a sequence of finite subsets of T satisfying Conditions (i) and 
(ii). Suppose a : R+ — > is a nondecreasing continuous function such that a(2h) < K\a(K) 
for some constant K\ > 0. If there exist constants K2 > 0, 7 G (0, 1] and Eq > such that 
limt7(/i)(log(l//i)) (1+£o)/7 = and 

y*l( max max \X(t p ) - X(t'-,)\A < K 2 a(2- n V . (2.1) 

Then X has a continuous version, still denoted by X, such that for all e > 

sup teT sup d( t)<h \X(t) - X(s)\ 

hm — — , w = 0, a.s. (2.2) 

a{h){\og{l/h)) {l+£)h 

Proof Given any s,t G D = Ufc^=i with d(s,t) < 2~ n , let {r p (s),n < p < q} and 
{r p (t),n < p < q} be the two approximating chains to s and t given by Condition (ii). The 
triangle inequality and the fact r n (s) = r n (t) imply 

1 1 
\X(t)-X(s)\< \X(r p (t))-X(r p ^(t))\+ ^ \X(r p (s)) - X(r p ^(s))\ 

p=n + l p=n+l 

OO v / 

<2 ™ , ™ ax . A x ( t p) - x ( t p-i)\- 

It is helpful to note that, for each p > n + 1, the maximum in (|2.3f) is taken over at most k,q N p 
increments, where N p denotes the cardinality of D p . 
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For any integer n > 1, let 

A n = sup sup \X(t) — X(s)\. 

s,teDd(s,t)<2- n 

It follows from ()2.3j) . the elementary inequality \a + 6| 7 < |o| 7 + |6| 7 (7 G (0, 1] and a, b G 
and (IO that 



E(A 7 )<2 7 V E^ max max |X(r„) - X(r' , ) l 7 ^) < 2 7 i^o"^ 



■n\-y 



p=n+l 

Hence Markov's inequality gives that for all integers n > 1 and real numbers it > 

P(A n > a(2~ n ) u^j < 2 7 K 2 u~ 7 . (2.4) 
For any e G (0, £0) and r/ > 0, by taking u = n{\og 2 n )( 1+£ )/ 7 , we derive from (|2.4p that 

A n > ^(2-") (log 2")( 1+£ )/ 7 ) <-^, 



which is summable. The Borel-Cantelli lemma implies almost surely 

A n < 7 ] a(2- n ) (log 2 n )( 1+£ )/ 7 (2.5) 

for all n large enough. Therefore, X(t) is a.s. uniformly continuous in D = U£Li -^fc- Since L> 
is dense in T, it is standard to derive from (|2.5p that X has a continuous version (still denoted 
by X) such that for all e > 0, 

sup tgT sup d( st)<2 -n -X(s)| 
hm V-r= ; , w = 0, a.s. (2.6) 

n^oo fj(2" ri ) (log2 n )( 1+£ )/ 7 

By (|2.6|) . the properties of <r and a monotonicity argument we derive (|2.2|) . This finishes the 
proof of Theorem 12.11 □ 

From now on, we will not distinguish X from its continuous version. The following result 
follows directly from Theorem 12.11 and is often more convenient to use. 

Corollary 2.2 Let X = {X(t),t G T} be a real-valued random field on a compact metric space 
(T,d) and let {D n ,n > 1} be as in Theorem \2. 11 Assume that there exist constants 7 G (0, 1], 
5 > and K > such that 



E 



(max max \X(r n ) - X^)] 1 ) <K2^ n (2.7) 



for all integers n > 1. Then for all e > 0, 

T su Pd(s,t) 

/,"(>gil /,))•'• " 



sup teT sup d( t)<h |X(t) -X(s)| 
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Next we apply Theorem 12.11 to a random field X = {X(t),t G [0,1]^} which may be 
anisotropic. Corollary 12.31 below is applicable to anisotropic stable random fields including lin- 
ear fractional stable sheets considered in Ayache, Roueff and Xiao (2007, 2009), harmonizable 
fractional stable sheets [cf, Xiao (2006, 2008)], operator-scaling stable fields with stationary 
increments in Bierme and Lacaux (2009) and others infinitely divisible fields. 

Given a constant vector (Hi, . . . , J3jv) £ (0> 1] j we consider the metric p on M. N defined 

by: 

p(s,t) = max \si -U\ H i Ms, teR N . (2.9) 

For every n > 1, define 

Dn = { • • • , , 1 < < j , VI < j < A (2.10) 



Then 

oo 

#D n <2 Qn and (J 5 n = [0, 1]". 

n=l 

In the above and in the sequel, j^D denotes the cardinality of D and Q = Ylj=i jf- ^ is easy 

to see that the sequence {D n ,n > 1} satisfies Conditions (i) and (ii) on the compact metric 
space ([0,l] N ,p). 

Corollary 2.3 Let X = {X(t),t G [0, 1]^} be a real-valued random field and assume a : — > 
]R+ satisfies the conditions of Theorem \2.1\ Let {D n ,n > 1} be the sequence defined above. If 
there exist constants 7 G (0, 1] and X > suc/i i/iaf 

VeT max max \X(t„) - X(r' i)! 7 ^ < K o(2~ n V 
for all integers n > 1, i/ien /or any e > 0, 

|X(t)-X(s)| 

lim sup sup ■ — ■ n , = 0, a.s. (2.11) 

«=[<),!]" |-t|<fc a(Ef=i - ^) I log (Ef=i l*i " «il H 0| ( 

Proof Note that the conditions of Theorem 12.11 are satisfied. By (|2.2p we have that for all 

lim ^o^s, t) < h \X(t)- X(s)\ = Q) ag (2 
a(h)(log(l/h)) — 

For any s,i€ [0,1]^, there is an integer n such that 2"™" 1 < p(s,i) < 2~ n . Then (f2T2j) . 
together with the properties of a, implies that 

\X(t)-X(s)\ <a(2- n )n^ <K 1 a(p(s,t)) ( log(l/p(«, t)j) ^ , a.s. 
for all s,i G [0, 1]^ such that p(s,t) is small. Since e > is arbitrary, this proves (|2.1ip , □ 
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3 Maximal moment index 



It is clear that condition (|2.ip [or (|2.7p ] is essential for Theorem 12.11 In many cases such as 
when X = {X(t),t G T} is a Gaussian or stable random field, or has stationary increments, 
one can normalize the random variables X(r n ) — X(r4_i) so that (|2.7p is reduced to conditions 
on the maximal 7-moments of a sequence of "homogeneous" random variables. 

Motivated by this, we introduce the notion of maximal 7-moment index for a (not necessar- 
ily stationary) sequence of random variables which, in turn, provides some sufficient conditions 
for (HZD to hold. 

Definition 3.1 Let k > 1} be a sequence of random variables such that, for some positive 
constants 7 and Ky, E(|£/c| 7 ) = Ky for all k > 1. Let M n {^) = E(maxi<fc< n |£fc| 7 ). Then the 
maximal 7 -moment (upper) index of{£k,k > 1} is defined by 

logM n ( 7 ) 

T = hmsup — . (3.1) 

twoo log n 

If 7 = 1, we simply call 9\ the maximal moment index of {£,k,k > 1}. 
Some remarks are in order. 

• If < 7 < (3, then Jensen's inequality implies that 6* 7 < X 9 p. 

• Since E(|£i| 7 ) < M n ( 7 ) < E(££ =1 l&l 7 )' we clearly have (9 7 G [0, 1]. 

• The maximal 7-moment index carries some information about the dependence structure 
and distributional properties of {£fc,/c > 1}. Usually the choice of 7 is determined by 
the heaviness of the tail probabilities of {£,k,k > 1}. For a fixed 7, smaller value of 
9y indicates more dependence among > 1}. This can be seen through the two 
extreme examples of stationary processes {^k,k > 1}: If = ^, then 9^ = 0; while if 
6c (k > 1) are i.i.d. and satisfy the conditions (|3.9p and (|3.1ip below, then # 7 = j/a for 
7 G (0,a). Further evidence can be found in Samorodnitsky (2004) where it is shown 
that the partial maxima of long-range dependent stable processes grow slower than those 
of short-range dependent processes. 

The following consequence of Theorem 12 . 1 1 shows the usefulness of maximal 7-moment index 
in determining the uniform modulus of continuity. 

Corollary 3.2 Let X = {X(t),t G 1^} be a continuous random field with values in M. and let 
{D n ,n> 1} C [0, 1] N be defined by 

^ = {(f'-"'^) : l<*i<2M<i<iv}. (3.2) 
We assume there exist constants 7 G (0, 1] and H > such that for all s, t G [0, 1]^ 

E(\X(s)-X(t)\~<) <K\s-t\ H "<. (3.3) 
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Consider the normalized random variables 
X(r p ) - X(r') 



[E(\X(t p )-X(t> )|7) 



1 1/7 ' 



V t p e Dp, r'_ j G O p _i(r p ) and V p > 1, (3.4) 



and number them according to the order D\, D2\D\, . . . and denote the sequence by k > 1}. 
{Cfc; > 1} ^ as a maximal ^-moment index 9 := # 7 and H*y > NO, then for every e > 0, 

r sup te[01]N sup ls _ tl<h \X(t) - X(s)\ 

limsup ^ = 0, a.s. (3.5) 

fc— 0+ h H ~~~ £ 

Namely, X{t) is uniformly Holder continuous on [0, 1]^ 0/ a// orders < H — 



Proof For any e > 0, it follows from (|3.ip that 

M 2 ^( 7 )=E( max |^r)<2^ + ^) 

for all n large enough. Combining this with f|3.3[) we derive 

Ef max max |X(r„) -^(T^i)! 7 ^ < K2^ (fl ~f " £)7 " 
V^eA 1 <_ 1 eo n -i(T„) / 

for all integers n large enough. Hence (|3.5p follows from Corollary 12.21 □ 

As an example, we show that the following multiparameter version of Kolmogorov's conti- 
nuity theorem follows from Corollary I3.2L 

Corollary 3.3 [Kolmogorov's continuity theorem] Let X = {X(t),t G [0,1]^} be a real- 
valued stochastic process. Suppose there exist positive constants (3, K and 5 such that 



E(\X(s) -Xit)] 13 ) < K\s-t\ N+s , V s,t G [0,1]". (3.6) 

13- 



Then X has a version which is uniformly Holder conditions on [0, 1]" of all orders < 4. 



Proof For any integer n > 1, let D n be defined by (|3.2p and let k > 1} be the sequence of 
the normalized random variables in (13. 4p . By (13. 6p , we see that (13. 3p is satisfied with 7 = f3 and 
i? = • In order to verify the second condition in Corollary 13.21 we note that E(|^| 7 ) = 1 
for every k > 1 and use the trivial bound E(max 1 </ C < m |6c| 7 ) — m for all integers m > 1 to 
derive # 7 < 1. Then it is clear that the conclusion of Corollary 13.31 follows from Corollary 13.21 
□ 

In the rest of this section, we study the maximal moment indices for several classes of 
random variables. The following lemma on Gaussian sequences will be useful in the next 
section. 

Lemma 3.4 Let {£fc,/c > 1} be a sequence of jointly Gaussian random variables with mean 
and variance 1. Then the following statements hold: 
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(i) There is a universal constant K± > such that 



E max |&| <K 4V %gii. (3.7) 

V l<fe<n / 

In particular, for every 7 £ (0, 1], the maximal j-moment index of k > 1} is 0. 

(n) Zf |E(^j^)| < 5 for a constant 5 £ (0, 1) and for all 1 < j < k < n, then there exists a 
constant such that E(max 1 </ C < n |^|) > K§ yTogn. 

Proof This lemma is well known, and we include a proof for completeness. Part (i) can be 
proved by using the metric entropy method. For any fixed integer n > 1, let T = {1, . . . ,n} 

1 /2 

equipped with the canonical metric d(i,j) = [E(£j — Cj) 2 ] ■ Then the d-diameter of T is 
at most 2. For any e G (0,1), the e-covering number Nd(T,e) < n. Hence Dudley's entropy 
theorem [cf. e.g., Marcus and Rosen (2006, Theorem 6.1.2) or Adler and Taylor (2007, Theorem 
1.3.3)] gives 

e( max \i k \] <K ^logn de, (3.8) 

V l<fe<71 / Jq 

which yields (|3.7p . 

Under the condition of (ii), we have d(i,j) = 2(1 — E(£j£j)) > 2(1 — 5) for all i 7^ j. Hence 
the conclusion of (ii) follows from the Sudakov minoration [see Talagrand (2006, Lemma 2.1.2)]. 
□ 

Dudley's entropy theorem has been extended to stochastic processes in Orlicz spaces, see 
Talagrand (2006, p. 30) for related references. Let : R+ — > R+ be a convex function such 
that ^(0) = and \&(r) > if r 7^ 0. For a random variable £, the Orlicz norm of £ is defined 

as 



q, = inf |c> : E*(^) < lj . 



If {£fc)^ > 1} is a sequence of random variables such that ||£fc||# = K < 00, then (1.53) in 
Talagrand (2006) implies that E( maxi<fc< n \£k\) — K^^in). This result extends Part (i) of 
Lemma [3. 4i By taking \I/(r) = r p (p > 1), one gets an upper bound for E(maxi<fc< n \£k\) when 
{£fcj k > 1} is a sequence of random variables with the same finite p-th. absolute moments. 
The following lemma is applicable to a-stable random variables. 

Lemma 3.5 Let {£fc,fc > 1} be a sequence of random variables. The following statements 
hold: 

(i) If there exist positive constants a and Kq such that 

P(|£fc| > u) < K 6 u~ a , V k > 1 and u > 0, (3.9) 
then for any 7 £ (0, a) and e > there is a finite constant Kj such that 

e( max \C k \A < K 7 n^ a (log n)( 1+£ ^/ a (3.10) 

\ l<fe<n / 

for all integers n > 2. Consequently, for any 7 < a, we have # 7 < j/ct. 



8 



(ii) If there exists a positive constant Kg such that 

> u) > K 8 u- a , Vfc>l and u>0 (3.11) 

and for all integers n > 1 and u > 



max |&| <u) < TT P(|&| <«), (3.12) 



fc=l 



i/ien, /or any 7 G (0, a), there is a constant Kg > such that 



E( max |£ fc n > Kgn^ (3.13) 



/or all integers n > 1. 



Remark 3.6 Clearly (|3.12p is satisfied if the random variables (/c > 1) are independent or, 
more generally, if the random variables \^ k \ (k > 1) are negatively orthand dependent, that is, 
for all n > 1 and x\, . . . , x n > 0, 



'(|£l| <*!,..., |£ n | <Xn) < Yl^{\^k\<X k ). 



k=l 



Moreover, by modifying the proof of Part (ii) of Lemma 13. 51 one can show that, if we assume 
in (i) that the random variables \^ k \ (k > 1) satisfies 



max |f fc | < u) > [JP(|^| <u), 



k=l 



then (I3.10P can be improved to E(maxi<fc<„ |£/c| 7 ) < Kn 1 ^ 01 . 

Proof of Lemma 13.51 Let 7 G (0, a) and e > be given constants. To prove Part (i), let 
$ : 1BL|_ — > be a nondecreasing and convex function satisfying <3?(0) =0 and 

x a/~f 

~ m— as x — > 00. (3.14) 

v y (logx) 1+e v y 

Here and in the sequel, f(x) ~ g(x) means f(x)/g(x) — ► 1 as x — > 00 or x — > 0. Then the 
inverse function of $(x), denoted by <J> _1 (x), is nonnegative, nondecreasing and concave on 
[0,oo). Moreover, 

$- l (x) ~ ^x 7 / Q (logx)( 1+E ) 7 / Q as x -> 00. (3.15) 
a 

By (J3I2D and (I3TT4"]) we derive E($(|£ fc |T)) < K w for all fe > 1, where K w is a positive 
constant depending on Kq, a, e and 7 only. This and Jensen's inequality together imply 



E( max \e k \A < 



< 



E$ max |£ fc P 



£ E (*(i&i 7 )) 

Kk<n 



<$>-\K w n) 
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Combining this and (|3.15|) yields (|3.10p . 
To prove Part (ii), we write 



E( max l^l 7 ) =7 f u 7 ~ x p( max \£ k \ > u) du. (3.16) 
It follows from (f3TT2]) and ([37TT]) that 



71 

'( max |£ fc | > u) > 1 - TT f 1 - P(|&| > 



fc=l (3.17) 



> 1 - ( 1 - K s u~ a 



for all u > Kl /a . By using the elementary inequality 



l-(l-x) n >-nx V < x < 1 - (i)V(n-i) 
we derive from f)3. lTj) that 

p( max > > -mi" Q (3.18) 

\ l<k<n ) 2 

for all u > K (l - (^) 1/(n_1) )" 1/a x n 1 /". Combining ([3TTB]) . (I3TT71) and dHlgjl we obtain 

E( max |^ fc | 7 J > / u 7 """ 1 ^ = Kn^ /a . (3.19) 
Vl<fc<n y y n i/ a 

This finishes the proof of Lemma 13.51 □ 

Part (ii) of Lemma l3.5l indicates that, in general, it is difficult to determine an optimal lower 
bound for E(maxi<fc< n |£fc| 7 ) without additional information about the dependence structure 
of {£fc,/c > 1}. We point out that, even for stationary stable sequences, it is an open prob- 
lem in general to determine sharp (i.e., up to constant factors) upper and lower bounds for 
E(maxi< fc < n |Cfc| 7 ). 

Here are some partial results. For symmetric stable or infinitely divisible sequences, some 
lower bounds can be derived from Theorem 2.3.1 and Theorem 5.3.2 in Talagrand (2006). For 
example, if is a stationary sequence of symmetric a-stable (SaS) random variables such 
that a £ (1,2] and there exists a constant 5 > such that ||£fc — Cj\\a > 8 for all k ^ j 
(where \\£\\ a denotes the scale parameter of £), then E( maxi<fe< n |^|) > K(\ogn)^ a ~ l ^/ a . 
On the other hand, Samorodnitsky (2004) studied the rate of growth of the partial maxima 
M n = max{|£fc| : 1 < k < n} of a stationary a-stable sequence {£,k,k > 1} based on the 
ergodic theoretical properties of the underlying flow. He discovered that (i) if the stationary 
SaS process {£fc, k > 1} is generated by a dissipative flow then M n grows always at the rate of 
n l l a and (ii) if the stationary SaS process k > 1} is generated by a conservative flow then 
M n grows at the rate slower than n}l a . Samorodnitsky (2004, p. 1440) conjectured that many 
other important properties of k > 1} will also change as the underlying flow changes from 
being dissipative to being conservative. In particular, we believe that, if X = {X(t),t £ M} is a 
(self-similar) SaS process with stationary increments, then the uniform modulus of continuity 
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of X depends on the nature of the flow generating the stationary sequence = X(k+1) — X(k) 
(k > 0). This idea will be pursued further elsewhere. 

Combining Lemma 13.51 with the proof of Theorem 12.11 we have the following result on 
modulus of continuity for general self-similar processes with stationary increments, which is 
an improvement of Corollary 13. 21 In the special case of N = 1 and X = {X(t),t G R} being 
an a-stable process, it recovers Theorem 2 in Kono and Maejima (1991b). 



Corollary 3.7 Let X = {X(t),t G R } be a real-valued H- self- similar random field with 
stationary increments. Let V = {v£, 1 < £ < 2 N — 1} be the set of vertices of [0, 1]^, excluding 
0. // there exists a constant a > jj such that 

F(\X(vi)\ > u) < Ku~ a , Vv e G V and u > 0, (3.20) 

then for any e > 0, 

sup ie[0 ,i]iv su P | s _ t | <A \X(t) -X(s)\ 

hmsup = = 0, a.s. (3.21) 

h H ~~ (logl//i)~ 



Proof Again we use the £°° metric in R^. For every n > 1, let D n be defined as in (|3.2p . 
Then the sequence {D n ,n > 1} satisfies the conditions in Section 2. The self-similarity and 
stationarity of the increments of X imply that, for every r n G D n and r^_i G O n _i(r ra ), there 
is a V£ G V such that 

X(r n ) - = |r n - r^il^X^). 

Hence we apply Lemma 13.51 to the normalized sequence { — -p^- ^^"h 1 ^ } to derive that for any 
< 7 < min{ 1 , a} 

e( max max |X(r„) - X^^lA < K 2^ H -^ n ( log 2^) ( 1+e ^ Q . (3.22) 
Hence (l3T2T]l follows from (I3T221) and Theorem [2TJ □ 



4 Applications to stable random fields 

In this section we apply the results in Sections [2] and [3] to harmonizable-type a-stable random 
fields with < a < 2. Similar methods can be applied to other types of stable random fields, 
or more generally, infinitely divisible processes. 

4.1 Harmonizable-type stable random fields with stationary increments 

Let X = {X(t), t G R^} be a real- valued stable random field defined by 

X(t) = Re / (e i{t ' x) - 1) M a (dx), (4.1) 
Jr n 
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where M a is a rotationally invariant a-stable random measure on M. N with control measure A, 
which satisfies 

/ (1 A \x\ a ) A(dx) < oo. (4.2) 

This condition assures that stochastic integral in (|4.1|) is well-defined, see Samorodnitsky and 
Taqqu (1994, Chapter 6) for further information. The measure A is called the spectral measure 
of X and its density function, when it exists, is called the spectral density of X. 

It can be verified that the stable random field X defined by (|4.ip has stationary increments 
and X(0) = 0. Denote the scale parameter of X(t) by Then for all t G R , 

\\X(t)\\% = 2 a l 2 f (l-cos(t,x)) a/2 A{dx). (4.3) 



Similar to Gaussian processes, this function plays an important role in studying sample path 
properties of stable random field X defined by (|4.ip . 

For simplicity, we will assume that the spectral measure A is absolutely continuous and its 
density function f(x) satisfies the following condition 

f(x) < K n \x\~ (aH+N) , Vx G M. N with |x| > K\2, (4.4) 

where Kn,K\2 > and H G (0,1) are constants. As shown by Theorem 14.51 below, the 
parameter H determines the smoothness of the sample function X(t). 

Now we provide some examples of stable random fields satisfying the condition (14.41) . 



Example 4.1 [Harmonizable fractional stable motion] Let H G (0, 1) and a G (0,2] be given 
constants. The harmonizable fractional stable field Z H = {Z H (t),t G R^} with values in R is 
defined by (|4.ip with spectral density 

(4.5) 



x 



where c(a, H, N) > is a normalizing constant such that the scale parameter of Z H (e\) equals 
1, where e x = (1, 0, . . . , 0) G R N . Hence we have \\Z H (t)\\ s = \t\ H for all t G R N . 

It is easy to verify that the a-stable random field Z H is .ff-self-similar with stationary 
and isotropic increments [or stationary increments in the strong sense in terms of Samorod- 
nitsky and Taqqu (1994, p. 392)]. The random field Z H is a stable analogue of fractional 
Brownian motion B H of index H and serves as an important representative for understanding 
harmonizable-type stable random fields. Even though the results on local times of Z H [Nolan 
(1989) and Xiao (2008)] and uniform modulus of continuity [see Theorem 14. 51 below] show that 
Z H shares many sample path properties with B H , little has been known about the sharpness of 
these results and the existing tools do not seem to be capable for attacking these problems. It 
is an interesting (and challenging) task to develop new methods for studying the fine structures 
of Z H and other stable random fields. 

Example 4.2 [Fractional Riesz-Bessel a-stable motion] Consider the stable random field 
X = {X(t),t G R^} in R defined by ((4~Tj) with spectral density 

J~(,V\ ) | x |2 7 ( 1 + | x |2)r,' V ) 
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where 77 and 7 are positive constants satisfying 



r? + 7>y, 0<2j <a + N 

and 0(0,7,77, N) > is a normalizing constant. When a = 2 [i.e., the Gaussian case] such 
density functions were consider by Anh et al. (1999). Since /~ „ involves both the Fourier 
transforms of the Riesz kernel and the Bessel kernel, Anh et al. (1999) called the corresponding 
Gaussian random field X the fractional Riesz-Bessel motion with indices 77 and 7. They 
showed that these Gaussian random fields can be used for modeling simultaneously long range 
dependence and intermittency. 

In analogy to the terminology in Anh et al. (1999), we call X the fractional Riesz-Bessel 
a-stable motion with indices 77 and 7. Clearly, the spectral density f 7>r) (x) in (|4.6p satisfies 
(|4.4p with H = (2(7/ + 7) — N)/a. Moreover, since the spectral density / 7i „(aj) is regularly 
varying at infinity of order 2(77 + 7) > N, by modifying the proof of Theorem 1 in Pitman 
|17j we can show that, if 2(7 + 77) — N < a, then ||X(t;)|| q is regularly varying at of order 
(2(77 + 7) - N)/a and 

||*(*)||a Ht| (2( ^ )-JV)/a > as |t|-»0. (4.7) 

We will see that the modulus of continuity of X is determined by (14. 7p . 

The following result provides information on the maximal moment of harmonizable-type 
stable random fields. It is more precise than Part (i) of Lemma 13.51 

Proposition 4.3 Let X = {X(t),t <E [0,1]^} be an a-stable random field defined by 1\ ) with 
spectral density satisfying \4-4\ )- Then for every < 7 < min{l, a}, r\ > and n > 1, 

e( max max \X(r n ) - Xir'^lA < K2^ n , (4.8) 

where {D n ,n > 1} is the sequence defined in \3.2i) . 

In order to prove Proposition 14.3^ we will make use of a LePage-type representation for X, 
which allows us to view X as a mixture of Gaussian process. This powerful idea was due to 
Marcus and Pisier (1984) and has become a standard tool for studying sample path regularity of 
stable processes. For more information on series representations of infinitely divisible processes 
and their applications, see Rosihski (1989, 1990), Kono and Maejima (1991a), Samorodnitsky 
and Taqqu (1994), Marcus and Rosinski (2005), Bierme and Lacaux (2009), just to mention a 
few. 

We need some notation. Let [i be an arbitrary probability on which is equivalent to 
the Lebesgue measure Denote by (p = dfi/dXjy the Radon-Nikodym derivative of /i, so 
fi(dx) = <p(x) dx. Assume that 

• j > 1} is a sequence of Poisson arrival times with intensity 1; 

• {gj, j > 1} is a sequence of i.i.d. complex valued Gaussian random variables such that 
g - £ e ie gj for all 6 € M and E(|Re#i| a ) = 1; 

• j > 1} is a sequence of i.i.d. random variables with values in and density function 
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• the sequences {Tj,j > 1}, {gj,j > 1} and > 1} are independent. We denote the 

expectations with respect to {Tj,j > 1}, {gj,j > 1} and > 1} by Ep, E 5 and E^, 

respectively. 

The following lemma is from Bierme and Lacaux (2009); see also Kono and Maejima (1991a) 
and Marcus and Pisier (1984). 

Lemma 4.4 For any family of complex-valued functions h(t, •) G ^"(R^, dx) (0 < a < 2), let 
Z = {Z(t),t £ M. N } be the a-stable random field defined by 



Z(t)=Re f h(t,x)Z a (dx), Vi G 



where Z a is a complex-valued, rotationally invariant a-stable random measure on (M. N , B(M )) 
with Lebesgue control measure. Then 

z(t),te R N \ = y(t),t g K 

where = means equality in finite dimensional distributions and 



Y(t) = CaReff^TJ^^y^h^&gX 

^ .7=1 ' 



(4.9) 



In the above, for every t G M. N , the random series converges almost surely and C a is the 

constant given by 

Proof of Proposition 14.31 Choose two positive constants n and (3 such that N — 2a ^S^ < 
(3 < N. Let tp : M^O} -> [0, oo) be the function defined by 

C ^ J #13 M"* if |*|<3, 

[ A 14 (|x| (log |x|) ^ ') it \x\ > 3, 

where the constants Ki% and i^i4 are chosen such that J^ N <p(x) dx = 1. 

By Lemma |4.4| the stable random field X defined by (|4.ip with spectral density /(a;) has 
the same finite dimensional distributions as 

Y(t) = C a Re (^rT 1 /^.)-!^^^.)^^ (4 U) 

where the function h(t,x) is defined by h(t,x) = {e 1 ^^ — l) f{x) l / a . Hence it is sufficient to 
prove (|4T8|) for Y. 
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Conditional on j > 1}, Y is a Gaussian random field with incremental variance 

given by 



E, 



(Y(t) - Y(s)) 2 ] =2C 2 a J2 r f a ^j)~ 2/a I 1 " cos <* " /fe) 2/C 



(4.12) 



where the inequality follows from (|4.4p and X > is a constant. 
Hence by using (|4.12p and Lemma 13.41 we have 



+-) 

a 1 



E, 



max max P^C^n) — ^( r n-i)| 7 ) 
T n eD„ ^eOfl-^Tn) / 



< IT max max 



¥. g {Y{r n )-Y{r' n _ x )f 



7/2 



7 /2 



(4.13) 



< if 



7 /2 



7/2 



It remains to show that (|4.8j) follows from taking expectations on both sides of (14.13P with 
respect to > 1} and {Tj,j > 1}. 



Note that, for every j > 1, Ij is a Gamma random variable with density function 



-1 „-x 



X? 6 

P(x) = _ , V x > 0. 

It is elementary to verify that there is a constant if > such that 

E r (iy 2/a ) <Kj- 2 / a , Vj>2/a 

and 



(4.14) 



E r (r\ 7/a ) < K V 1 < j < 2/a. (4.15) 
Since 7 £ (0, 1), we use the elementary inequality (x + y) 7 ^ 2 < £ 7 / 2 + y 7 ^ 2 and Jensen's 



inequality to derive 

E r,e| (E r 7 2/ Vfe)" 2/a {1 a 2- 2 %f} 

[2/a 



T /2' 



/ L Z /«J \ 

<E r)5 ( £ rT 7 /vfe)- 7/Q {iA2- 7 "ic J r}ie,r 7( ^ ) 

V 3=1 ' 



(4.16) 



+ 

:=/i+/ 2 - 



Er, ( £ r 7 2/ V(^)" 2/Q {1 A 2"*%| 2 } |Ci|- 2(if+ ^ 

j=L2/«J + l 



17/2 
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Since I\ and I2 can be estimated by using the same method, we only consider I2 below. 
Taking the expectation Eg first, we have 



e € (^)- 2/q {ia2- 2 "|^i 2 }i^i- 2(h+ ^ ) ) 

^Oe) 1 -^ (1 A 2~ 2 >| 2 ) 



dx 



\x\ 2 ( H+ ^ 



(4.17) 



2 21 

\x\<2 n 
Jl + J 2 - 



dx 



.|2(ff+£-l) J ]xl>2 



(fix) 1 a 



dx 



\x\ 2 ^ H+ ^ 



By ()4.10p and a change of variables, we derive 



Ji = K2 



~2n 



dr 



Q r 2H+{N-f5){l-l)-l 

2 " dr 



+ 



3 r 

2 



2J/-1 



log ry- 1 



(4.18) 



< K2- 2 ^" n-( 1+ ")( 1 -^). 
Note that, because of the choice of /?, the first integral is constant. Similarly, we also have 



j K2 -2Hn n -(l+r,)(l-i)_ 



(4.19) 



By (liTH) . (jUBD and gJJJ]), we obtain 



/ 2 < 



j=[2/a\+l 

< ET2-^ n n- 7(1+r ' )( 5-i). 
Similarly, we can derive 



/ 00 \ -1 7/2 

E t (k T- 2/a 2- 2Hn n^ 1+ ^ 1 -^\ 



j=[2/a\+l 



-2/a 



7/2 



(4.20) 



Combining ([4TTU]) . (|4T2"Uj) and (pOT]) we obtain 

E r, 5 | (Er7 2/ Vfe)- 2/Q {1 A 2- 2 %f} \^ H+ ^ 

Finally (fl~8|) follows from (f4"7T3|) and (ET22D . This proves Proposition S3J 



(4.21) 



7/2 ' 



(4.22) 



□ 
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The following is a consequence of Proposition 14.31 and Theorem 12.11 

Theorem 4.5 Let X = {X(t),t £ R^} be an a-stable random field defined by fa4-l\ ) with 
spectral density satisfying \4-4\ )- Then for any e > 0, 

sup te[01]N sup ls _ tl<h \X(t)-X(s)\ 

hmsup ■ — - ■ t— . . = 0, a.s. (4.23) 

h H {logl/h) (2+£)/a 

Proof It follows from Proposition 14.31 that for all rj > and integers n > 1, 



VEf max max \X(t v ) - X(t' ,)\A < K 2- H ^ n 



7(1+3) 



Thus, X satisfies (JXT]) with cx(/i) = /^(log 1/h)— . Since ry > and < 7 < a are arbitrary, 
(|4T23|) follows from Theorem □ 



4.2 Harmonizable fractional stable sheets 

For any given < a < 2 and H = (Hi, ...,Hn) G (0,1) , we define the harmonizable 
fractional stable sheet Z H = {Z H (t),t 6 R^} with values in R by 

Z & {t) = Re / J] ^ ~ r Z a (dA), (4.24) 

where Z a is a complex- valued random measure as in Lemma 14.41 

From (|4.24p it follows that Z H has the following operator-scaling property: For any N x N 
diagonal matrix E = (b(j) with bu = b% > for all 1 < i < N and bij = if i 7^ j, we have 

{Z ff (Et),teR N } = ^(jJbf j ^jZ a (t), teR^j. (4.25) 

Along each direction of R^X, Z H becomes a real- valued harmonizable fractional stable motion 
[cf. Samorodnitsky and Taqqu (1994, Chapter 7)]. When the indices H\, . . . , -Hjv are n °t the 
same, Z H has different scaling behavior along different directions and this anisotropic nature 
induces some interesting geometric and analytic properties for Z . Note that Z H does not 
have stationary increments in the ordinary sense and, thus, is different from the operator- 
scaling stable fields considered in Bierme and Lacaux (2009). See Xiao (2006, 2008) for further 
information on sample path properties of stable random fields. 

The following result gives the uniform modulus of continuity for Z H , which is significantly 
different from the result for linear fractional stable sheets obtained in Ayache, Roueff and Xiao 
(2007, 2009). 

Theorem 4.6 For any arbitrarily small e > 0, one has 

\Z**(s) — Z^(t)\ 

lim sup — 1 ^ — — tt— j- = 0, a.s. (4.26) 

^*6|P,1F,|-*|<A Ef=i \*i ~ tj\ H > I log ( Ef=i \sj ~ |(2+E)/Q 



"3\ \ 
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Proof Let p be the metric on M> N defined in (|2.9p and let {D n ,n > 1} be the sequence 
defined in (|2.10p . We claim that for all integers n > 1, < 7 < min{l, a} and 77 > 0, 

e( max max |Z*(t„) - ^ 5 «_ 1 )| 7> ) < i^^ra 1 ^. (4.27) 

Before proving (|4.27p . we note that it implies 

00 , _ _ \ 

VE max max |Z^(r p ) - Z S (t' ,) I 7 ) < K 2^ n n 11 ^. 

^ \r p €D p r;-ieOp-i(r p ) ' / 

Hence satisfies the conditions of Corollary O with a(h) = h{ log l/ti) {1+v)/a . Thus, (1061) 
follows from (12, lip . 

It remains to prove (|4.27p . Since this is similar to the proof of Proposition 14,31 we only 
provide a sketch of it. Let 77 and (3 be two positive constants with (3 satisfying 



/ 2a(l - Hj) \ 

max < 1 — > 

i<j<N { 2 - a J 



< /3 < 1. 



Define the density function tp : M. N — > [0, 00) by ^(a?) = Y\f=i <Pj( x j)) where x = (x±, . . . , xjy) 
and 

J i^i 5 |xj| _/3 if < 3, 

Wi) ~ \ K W (^KlogN)^)- 1 if \ Xj \ > 3. 

In the above, the constants and K\§ are chosen such that J" RJV tp(x) dx = 1. By Lemma 
Ol has the same finite dimensional distributions as 

Y(t) = C a Jte(f^rj 1 ' a <p{Z j )-V«h(t,Z j )g\ 

v j=1 / 

where the function h(t,x) is now defined by 

Again, by conditional on {(£,,1^), j > 1}, we have 



E 9 



-y(,)) 2 ] = c 2 XX 2A V&r 2/Q |*(t,&)-fc(,,&)| a . (4.28) 



If follows from (|4.28[) and Lemma 13.41 that 
El max max I^C^n) — ^( r n-i) | 7 ) 

< if n 7 / 2 Erf| f fr^Vfe)- 2 /" max max 



|^(Tn,^') - k( r n-l>£j) 



7 /2- 



(4.29) 
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It can be verified that for every s, t and x £ 



\h(t, x) — h(s, x)\ 



N N 2 N 

TT ( e « fc *i, - 1) - TT (e iSfea:fe - 1) T7 



fc=i 



fc=i 



fc=l Ffc 



A? 



< 



{lA(t e -s e ) 2 x 2 } H{lAt 2 k x 2 k } H{lAs 2 k xl} 



k>£ 



N 1 



(4.30) 



k=l \ x k 



By using (14.30f) and an argument similar to the proof of Proposition 14.31 one can derive 



E rA E r 7 2A V&)- 2/a max |^r n , - 



7/2 ' 



< K2" 7n n" 7(1+, ' )( 5-^). 
This, together with (|CZ5|> . proves (|Q7)) . 



□ 



References 

[1] Adler, R. J. and Taylor, J. E. (2007), Random Fields and Geometry. Springer, New York. 

[2] V. V. Anli, J. M. Angulo and M. D. Ruiz-Medina (1999), Possible long-range dependence in 
fractional random fields. J. Statist. Plann. Inference 80, 95-110. 

[3] Ayache, A., Roueff, F. and Xiao, Y. (2007), Local and asymptotic properties of linear fractional 
stable sheets. C. R. Acad. Sci. Paris, Ser. A. 344, 389-394. 

[4] Ayache, A., Roueff, F. and Xiao, Y. (2009), Linear fractional stable sheets: wavelet expansion 
and sample path properties. Stoch. Process. Appl. 119, 1168-1197. 

[5] Bernard, P. (1970), Quelques proprietes des trajectoires des fonctions aleatoires stables sur R k . 
Ann. Inst. H. Poincare Sect. B (N.S.) 6, 131-151. 

[6] Bierme, H. and Lacaux, C. (2009), Holder regularity for operator scaling stable random fields. 
Stoch. Process. Appl. to appear. 

[7] Csaki, E. and Csorgo, M. (1992), Inequalities for increments of stochastic processes and moduli 
of continuity. Ann. Probab. 20, 1031-1052. 

[8] Garsia, A. M. (1972), Continuity properties of Gaussian processes with multidimensional time 
parameter. Proc. 6th Berkeley Symp. Math. Statist. Probability II, pp. 369-374, University of 
California Press, Berkeley. 

[9] Khoshnevisan, D. (2002), Multiparameter Processes: An Introduction to Random Fields. 
Springer, New York. 

[10] Kono, N. and Maejima, M. (1991a), Self-similar stable processes with stationary increments. In 
Stable processes and related topics (Ithaca, NY, 1990), 275-295, Progr. Probab., 25, Birkhauser 
Boston, Boston, MA. 

[11] Kono, N. and Maejima, M. (1991b), Holder continuity of sample paths of some self-similar stable 
processes. Tokyo J. Math. 14, 93-100. 



19 



[12] Kwapieh, S. and Rosinski, J. (2004), Sample Holder continuity of stochastic processes and ma- 
jorizing measures. In: Seminar on Stochastic Analysis, Random Fields and Applications IV, pp. 
155-163, Progr. Probab., 58, Birkhauser, Basel. 

[13] Marcus, M. B. and Pisier, G. (1984), Characterizations of almost surely continuous p-stable 
random Fourier series and strongly stationary processes. Acta Math. 152, 245-301. 

[14] Marcus, M. B. and Rosen, J. (2006), Markov Processes, Gaussian Processes, and Local Times. 
Cambridge University Press, Cambridge. 

[15] Marcus, M. B. and Rosinski, J. (2005), Continuity and boundedness of infinitely divisible pro- 
cesses: a Poisson point process approach. J. Theoret. Probab. 18, 109-160. 

[16] Nolan, J. (1989), Local nondeterminism and local times for stable processes. Probab. Th. Rel. 
Fields 82, 387-410. 

[17] Pitman, E. J. G. (1968), On the behavior of the characteristic function of a probability sidtri- 
bution in the neighbourhood of the origin. J. Australian Math. Soc. Series A 8, 422-443. 

[18] Rosinski, J. (1989), On path properties of certain infinitely divisible processes. Stochastic Proc. 
Appl. 33 , 73-87. 

[19] Rosinski, J. (1990), On scries representations of infinitely divisible random vectors. Ann. Probab. 
18, 405-430. 

[20] Samorodnitsky, G. (2004), Extreme value theory, ergodic theory and the boundary between short 
memory and long memory for stationary stable processes. Ann. Probab. 32, 1438-1468. 

[21] Samorodnitsky, G. and Taqqu, M. S. (1994), Stable non-Gaussian Random Processes: Stochastic 
m,odels with infinite variance. Chapman & Hall, New York. 

[22] Takashima, K. (1989), Sample path properties of ergodic self-similar processes. Osaka Math. J. 
26, 159-189. 

[23] Talagrand, M. (2006), Generic Chaining. Springer- Verlag, New York. 

[24] Xiao, Y. (2006), Properties of local nondeterminism of Gaussian and stable random fields and 
their applications. Ann. Fac. Sci. Toulouse Math. XV, 157-193. 

[25] Xiao, Y. (2008), Properties of strong local nondeterminism and local times of stable random 
fields. Seminar on Stochastic Analysis, Random Fields and Applications VI, to appear. 



Yimin Xiao. Department of Statistics and Probability, A-413 Wells Hall, Michigan State 

University, East Lansing, MI 48824, U.S.A. 

E-mail: xiao@stt.msu.edu 

URL: http : //www . stt . msu . edu/~xiaoyimi 



20 



